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Abstract 


We show that initial-state interactions contribute to the cos 2</> distribu¬ 
tion in unpolarized Drell-Yan lepton pair production pp and pp —> 
without suppression. The asymmetry is expressed as a product of chiral-odd 
distributions hi{xi,p\) x hi {x 2 ,k‘]_), where the quark-transversity function 
hi{x,p‘]_) is the transverse momentum dependent, light-cone momentum dis¬ 
tribution of transversely polarized quarks in an unpolarized proton. We com¬ 
pute this (naive) T-odd and chiral-odd distribution function and the resulting 
cos 2(j) asymmetry explicitly in a quark-scalar diquark model for the proton with 
initial-state gluon interaction. In this model the function hi{x,p\) equals the 
T-odd (chiral-even) Sivers effect function f^{x,p\). This suggests that the 
single-spin asymmetries in the SIDIS and the Drell-Yan process are closely re¬ 
lated to the cos 2(^ asymmetry of the unpolarized Drell-Yan process, since all 
can arise from the same underlying mechanism. This provides new insight re¬ 
garding the role of quark and gluon orbital angular momentum as well as that 
of initial- and final-state gluon exchange interactions in hard QCD processes. 


1 Introduction 

Single-spin asymmetries in hadronic reactions have been among the most challenging 
phenomena to understand from basic principles in QCD. Several such asymmetries 
have been observed in experiment, and a number of theoretical mechanisms have 
been suggested 1|, 0 B i B §. Recently, a new way of producing single-spin asym¬ 
metries in semi-inclusive deep inelastic scattering (SIDIS) and the Drell-Yan process 
has been put forward |^, ||. It was shown that the exchange of a gluon, viewed as 
initial- or hnal-state interactions, could produce the necessary phase leading to a sin¬ 
gle transverse spin asymmetry. The main new feature is that despite the presence 
of an additional gluon, this asymmetry occurs without suppression by a large energy 
scale appearing in the process under consideration. It has been recognized since then 
[^, that this mechanism can be viewed as the so-called Sivers effect i0. which 
was thought to be forbidden by time-reversal invariance [Q. Apart from generating 
Sivers effect asymmetries, the mechanism offers new insight regarding the role of or¬ 
bital angular momentum of quarks in a hadron and their spin-orbit couplings; in fact, 
the same S ■ L matrix elements enter the anomalous magnetic moment of the proton 
[0. The new mechanism for single target-spin asymmetries in SIDIS necessarily re¬ 
quires non-collinear quarks and gluons, and in the Sivers asymmetry the quarks carry 
no polarization on average. As such it is very different from mechanisms involving 
transversity (often denoted by hi or (5g), which correlates the spin of the transversely 
polarized hadron with the transverse polarization of its quarks. 

In further contrast, the exchange of a gluon can also lead to transversity of quarks 
inside an unpolarized hadron. This chiral-odd partner of the Sivers effect has been dis¬ 
cussed in Refs. 10, and in this paper we will show explicitly how initial-state inter¬ 
actions generate this effect. Goldstein and Gamberg reported recently that hi{x,p‘]_) 
is proportional to f^{x,p\) in the quark-scalar diquark model [^. We confirm this 
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and find that these two distribution functions are in fact equal in this model. Al¬ 
though this property is not expected to be satished in general, nevertheless, one may 
expect these functions to be comparable in magnitude, since both functions can be 
generated by the same mechanism. We investigate the consequences of the present 
model result for the unpolarized Drell-Yan process. We obtain an expression for the 
cos 20 asymmetry in the lepton pair angular distribution. Here 0 is the angle between 
the lepton plane and the plane of the incident hadrons in the lepton pair center of 
mass. This asymmetry was measured a long time ago ilT^, ITil and was found to be 


large. Several theoretical explanations (some of which will be briefly discussed below) 
have been put forward, but we will show that a natural explanation can come from 
initial-state interactions which are unsuppressed by the invariant mass of the lepton 
pair. 


2 The unpolarized Drell-Yan process 


The unpolarized Drell-Yan process cross section has been measured in pion-nucleon 
scattering: 7i~ N —>■ fi~ X, with N deuterium or tungsten and a 7r“-beam with 
energy of 140, 194, 286 GeV 


131 and 252 GeV |14]. Gonventionally the differential 


cross section is written as 
1 da 


— = — - f 1 -I- A cos^ 6 + li sin^ 6 cos (p + — sin^ 6 cos 20^ . 

(TaiZ47rA-|-3\ 2 / 


( 1 ) 


These angular dependencies can all be generated by perturbative QGD corrections, 
where for instance initial quarks radiate off high energy gluons into the hnal state. 
Such a perturbative QGD calculation at next-to-leading order leads to A ~ 1,/r ~ 
0, z/ 0 at very small transverse momentum of the lepton pair. More generally, the 

pA[| is expected to hold at order and the relation 
perturbative QGD corrections 

T1 


Lam-Tung relation 1 — A — 2z^ = 0 


is hardly modihed by next-to-leading order 
However, this relation is not satished by the experimental data 
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The Drell-Yan 


data shows remarkably large values of z/, reaching values of about 30% at transverse 
momenta of the lepton pair between 2 and 3 GeV (for = (4 — 12 GeV)^ and 

to be discussed below). These large values 
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extracted in the Gollins-Soper frame 
of ly are not compatible with A ~ 1 as also seen in the data. 

A number of explanations have been put forward, such as a higher twist effect 
20| , ^ , following the ideas of Berger and Brodsky . In Ref. the higher twist 


effect is modeled using an asymptotic pion distribution amplitude, and it appears to 
fall short in explaining the large values of z/. 

In Ref. factorization-breaking correlations between the incoming quarks are 
assumed and modeled in order to account for the large cos 20 dependence. Here the 
correlations are both in the transverse momentum and the spin of the quarks. In 
Ref. H this idea was applied in a factorized approach involving the chiral-odd 


partner of the Sivers effect, which is the transverse momentum dependent distribution 


^We neglect sin0 and sin 20 dependencies, since these are of higher order in Ug 00 and 
expected to be small. 


are 
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function called h^. From this point of view, the large cos 20 azimuthal dependence can 
arise at leading order, i.e. it is unsuppressed, from a product of two such distribution 
functions. It offers a natural explanation for the large cos 20 azimuthal dependence, 
but at the same time also for the small cos0 dependence, since chiral-odd functions 
can only occur in pairs. The function is a quark helicity-flip matrix element and 
must therefore occur accompanied by another helicity flip. In the unpolarized Drell- 
Yan process this can only be a product of two functions. Since this implies a change 
by two units of angular momentum, it does not contribute to a cos 0 asymmetry. In 
the present paper we will discuss this scenario in terms of initial-state interactions, 
which can generate a nonzero function h^. 

We would also like to point out the experimental observation that the cos 20 
dependence as observed by the NAIO collaboration does not seem to show a strong 
dependence on A, i.e. there was no signihcant difference between the deuterium and 
tungsten targets. Hence, it is unlikely that the asymmetry originates from nuclear 
effects, and we shall assume it to be associated purely with hadronic effects. We refer 

m 


to Ref. 


for investigations of nuclear enhancements. 


We compute the function and the resulting cos 20 asymmetry explicitly 

in a quark-scalar diquark model for the proton with an initial-state gluon interac¬ 
tion. In this model hi{x,pY) equals the T-odd (chiral-even) Sivers effect function 
f^{x,pY). Hence, assuming the cos 20 asymmetry of the unpolarized Drell-Yan pro¬ 
cess does arise from nonzero, large this asymmetry is expected to be closely 
related to the single-spin asymmetries in the SIDIS and the Drell-Yan process, since 
each of these effects can arise from the same underlying mechanism. 

The Tevatron and RHIC should both be able to investigate azimuthal asymmetries 
such as the cos 20 dependence. Since polarized proton beams are available, RHIC 
will be able to measure single-spin asymmetries as well. Unfortunately, one might 
expect that the cos20 dependence m pp ^ ^^X (measurable at RHIC) is smaller 
than for the process 7r~ N ^ since in the former process there are no 

valence antiquarks present. In this sense, the cleanest extraction of would be from 
pp^ilX. 


3 Cross section calculation 

In this section we will assume nonzero and discuss the calculation of the leading 
order unpolarized Drell-Yan cross section (given in Ref. 0 with slightly different 
notation) 


da{hih2 —> iiX) 


dVLdxidx2d?q^ 




a,a 


B{y) cos(20) X 


(2h-pj_ h k^ — p^ kj 


h^hf 

M1M2 


( 2 ) 
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This is expressed in the so-called Collins-Soper frame ||19||, for which one chooses the 
following set of normalized vectors (for details see e.g. ||25||): 


i 

z 

h 


q/Q, 

q±/Q± = {q-XiPi -X 2 P 2 )/Q±, 


( 3 ) 

( 4 ) 

( 5 ) 


where Pi = Pi — q/{2xi), Pi are the momenta of the two incoming hadrons and q is 
the four momentum of the virtual photon or, equivalently, of the lepton pair. This 
can be related to standard Sudakov decompositions of these momenta 




2x1 


n^ + 


xiMl 


■ 


3^2 ^ Q 

2Q 


W 


2X2 


■ 


j ^ n>^ + q±, 


( 6 ) 

( 7 ) 

( 8 ) 


with Q\ = —qj_ = -c g^, via the identification of the light-like vectors 


= 


= 


P + F - 


( 9 ) 

( 10 ) 


The azimuthal angles lie inside the plane orthogonal to t and z. In particular, dO, 
= 2dydcl)\ where 0^ gives the orientation of = {g^'^ — FF -(- ly, the perpen¬ 

dicular part of the lepton momentum /; 0 is the angle between h (the direction of q^) 
and /_L. In the cross sections we also encounter the following functions oi y = l~/q~, 
which in the lepton center of mass frame equals y = {1 + cos 0)/2, where 6 is the angle 
of the momentum of the outgoing lepton I with respect to 5 (cf. Fig. |lD: 

A{y) = + + ’ ( 77 ) 

B{y) = y{l-y) = ^sm‘^e. ( 12 ) 

Furthermore, we use the convolution notation 


P 


ff 


j d^p^ d^k^ 5‘^{p^ + k^- qr^)/“(A,pl)/“(A, k\), 


( 13 ) 


where A, A are lightcone momentum fractions and a is the flavor index. 
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Figure 1: Kinematics of the Drell-Yan process in the lepton center of mass frame. 


In order to obtain the cross section expression one contracts the lepton tensor 


with the hadron tensor [g, ^ 
1 


+ 


q ^ —q 


= IJ dk^d‘^p^(fk^6^{p^+k^~q^) Tr (<F(p) 7 ^ $(fc) 7 ^) 

(14) 

where = AP^ = Aq^/xi, k = AP 2 = Aq 1x2- The correlation function $ is 
parameterized in terms of the transverse momentum dependent quark distribution 
functions in 


M V P Qo- p^_ 


(15) 


-hiT(A,r_L 




M 




M 2 


and similarly for 4). 

We end this section by giving the resulting expression for v 


= 2 E^a-^ 


(2h-pj_ h kj^ — Pj^ kj 


hi hi 

Ml M2 




(16) 


4 Asymmetry calculation 

The above cross section in terms of $ and <F can be represented by the diagram in 

Fig. 0. _ 

Insertion of the parameterization of <F and <F will yield the cos2(;/) asymmetry, 
among many other terms. However, in the lowest order quark-scalar diquark model 
the diagram Fig. ^ will not lead to nonzero hf- in <F, and consequently, also not 
to a nonzero cos 20 asymmetry. To generate such an asymmetry we will include 
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Figure 2: The leading order contribution to the Drell-Yan process 



Figure 3: The initial-state interaction contribution to the Drell-Yan process. 
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initial-state interactions corresponding to diagrams such as those depicted in Fig. 
^ Following the reasoning of Refs. this should be equivalent to Fig. g with 

an effective <F (and <h) with nonzero function. Here we do not intend to give a 

full demonstration of this in the Drell-Yan process; a generalized factorization theo¬ 
rem which includes transverse momentum dependent functions and initial/hnal-state 
interactions remains to be proven [^|. Instead we present how to arrive at an effec¬ 
tive <I> from initial/hnal-state interactions and use this effective $ in Fig. ^ Also, 
for simplicity we will perform the explicit calculation in QED. Our analysis can be 
generalized to the corresponding calculation in QCD. The hnal-state interaction from 
gluon exchange has the strength —>• where e* are the photon couplings 

to the quark and diquark. 

The diagram in Fig. coincides with Fig. 6(a) of Ref. used for the evalua¬ 
tion of a twist-4 contribution (~ 1/Q^) to the unpolarized Drell-Yan cross section. 
The differences compared to Ref. are that in the present case there is nonzero 


transverse momentum of the partons, and the assumption that the matrix elements 
are nonvanishing in case the gluon has vanishing light-cone momentum fraction (but 
nonzero transverse momentum). This results in an unsuppressed asymmetry which 
is a function of the transverse momentum Q± of the lepton pair with respect to the 
initial hadrons. If this transverse momentum is integrated over, then the unsuppresed 
asymmetry will average to zero and the diagrams will only contribute at order 1/Q^ 
as in Ref. |2^. 

First we will calculate the <I> matrix to lowest order (called 4*2^) in the quark- 
scalar diquark model which was used in Ref. [^. (Although the model is based on 
a point-like coupling of a scalar diquark to elementary fermions, it can be softened 
to simulate a hadronic bound state by differentiating the wavefunction formally with 
respect to a parameter such as the proton mass.) As indicated earlier, no nonzero 
and will arise from Next we will include an additional gluon exchange 
to model the initial/hnal-state interactions (relevant for timelike/spacelike processes) 
to calculate and do obtain nonzero values for and h/. Our results agree 
with those recently obtained in the same model by Goldstein and Gamberg []^. We 
can then obtain an expression for the cos 20 asymmetry from Eq. (|^) and perform 
a numerical estimation of the asymmetry. 


4.1 $ matrix in the lowest order 

As indicated in Fig. ^the initial proton has its momentum given by = {P^, P~, P±) 
(P+, ^,0^), and the hnal diquark P'^ = (P'+, P'", P'^) = (P+(l-A), ^ ^±)- 

We use the convention ± a^, a ■ b = |(a+6“ -|- a“6+) — a± ■ b±. 

We will hrst calculate the <I> matrix to lowest order ($2^) in the quark-scalar 
diquark model used in Ref. [0]. By calculation of Fig. H one readily obtains 


4>f = 


ag 


u{P,S)- 


f + m i/3r ■/ + m 


Lr 




MP, s) 


P+(l-A) 


(17) 
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Figure 4: 


Diagram which gives the lowest 


order $ (called $ 2 ^. 


= ag 


0^, 


m(P, S')(/+m) + m)u{P,S) 


P+(l - A) 


X 


^a(M2 - 


A 


with a constant a = l/(2(27r)^). The normalization is hxed by the condition 


J dA dVi /i(A, ri) = 1 . 


(18) 


In Eq. (P!7D we used the relation 


= r^(r' — 


= A{M^ 


r 

y 2 _|_ ^2 ^2 _|_ ^2 


) = r+{p-- 


X^ + rl 
(1 - A)P^ 


, 


1 - A 


A 


(19) 


This model is similar to the so-called spectator model (see e.g. Ref. [^), where 
in addition a vector diquark is included and the coupling constant g is treated as a 
form factor (in order to guarantee convergence). Of course, this can be assumed in 
the present model calculation as well and will be discussed in Section Assuming 
real form factors, the functions and are strictly zero in the spectator model. 


4.1.1 Calculation of/i( A, rj^) 

For the calculation of the denominator of the asymmetry one needs to know the 
function /i(A,r^), which can be obtained from given in Eq. (P!5|): 

A(A.rx) = iEit“'’(7Y“. (20) 

^ ±s ^ 


We now take <F = <I>l and for the numerator spinor contraction, we calculate 


1 

2 


^\u{P, S){f + m) ^ +m)u{P,S) ( 7 +) 


±s 


\/3q 
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( 21 ) 


= -Tr 
2 

= 2P+ 


(P + M) {f + m)'y~^{f+m) 
r]_ + (AM + mY 


Then, from Eqs. ( 0 ), (^) and (^If) , we arrive at 

/i(A,r_L) = + (AM + m)'"^ 


(- 


(1-A) ^a(M2 
(1 — A) r]_ + (AM + m)^ 


m2+r| 


A2+r 2 
l-A 


2(27r) 


+ By 


= C- 




+ D 


+ B) 


2 ’ 


where we define C = g‘^{l — A)/(2(27r)^), D = (AM + m)^ and 

A2 . 


B = A(l-A)(-M" + 


m 


+ 


A 1 -A' 


( 22 ) 


(23) 


Since we consider the proton state with mass M as a bonnd state composed of a 
qnark with mass m and a diquark with mass A, the function B as given in Eq. 
(|25|) is always nonzero and positive.The integral in Eq. (P!B[ ) with /i(A,r^) given in 
Eq. (P^) can for instance be regulated by assuming a cutoff in the invariant mass: 




k I i+m; 


^ < A^, and the value of is adjusted to satisfy the normalization 


condition Eq. 
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4.2 $ matrix with final-state interaction 

In order to obtain the $ matrix with final-state interaction (called from which 

one can trivially obtain the one with initial-state interaction, we calculate the diagram 
given in Fig. §(b). This is equal to the diagram calculated by Ji and Yuan to 
obtain nonzero /j^, starting from the formal gauge invariant definition of this trans¬ 
verse momentum dependent distribution function [P, In Fig. |](b) we attached the 
virtual photon line to the later end of the eikonal line in order to emphasize that the 
final-state interaction effect has become an ingredient of the distribution functions of 
the target proton. In reality, the whole eikonal line should be considered to be at the 
same point. 

Dehning through Fig. |^(b) (in the Feynman gauge), we have 




= lag 6162 


(24) 


d^k 1 -P+((l-x) + (l-A)) 
(27r)4 2 P+(l-A) 


X 


X 


P+((x — A) -I- ie) 




u{P,S){f+m) {^ + m)u{P,S) 


(/c^ —m? + ie){{k — r)^ — -A ie){{k — Py — A^ -|- ie) — w? 
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q 



- (h.c.) 

-iag^ 6162 
2 


(Pkx 

2(27r) 


P^dx 


(-(l-x)-(l-A)) 


P+3 X (x - A) (1 - x) 2P+(1 - A) 


X 


((x - A) + ie) a(M2 - 


X dk 


u{P, S){f + m) 


/^r 


+ m)u{P, S) 
1 


X 


k- - 
(h.c.) , 


I ( ( j .. 




-P-) + 


{\^+k^)-ie 

(l-x)P+ 


where we used /c+ = xP+. The derivation of the starting formula of Eq. (|2^) is given 
in the Appendix. This underlies the step from Fig. |^(a) to Fig. |^(b) and hence the 
step from Fig. ^ to Fig. 

For in Eq. (|^) , we consider only the contribution from the imaginary part 
of ((j;_A)+^e) ’ f^3at is, the contribution from —i 7 r 6 {x — A). There is no contribution 

from the real part of since the hermitian conjugate term cancels it. Then, 
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we have 


= 4 {-iag^ 6162 ) J 

X — in6{x —/S) 
X dk 


d?k± 

2(27r) 


J P+dx 


(25) 

-(l-x)-(l-A)) 


P+^ X (x-A) (l-x) 2P+(1-A) 


1 


A(M2 - A!±ri _ ^!±Il 


1-A 




X 


m(P, S')(/+m) (/(f + m)M(P, S') 

1 


7 _ _ (m2+fc^)-te 
xP+ 


r~) 


(A2+(fc_L-'rx)2)-ie 


(a:-A)P+ 

When we perform the k~ integration, we have 


p-\ I {X^+kj_)-i£ 
^ ^ il-x)P+ 


<s>f 


(26) 


= 4 TT 6162 


x6{x — A) 






((l-x) + (l-A)) 


2(27r)4y P+3 x (a: - A) (1 - a:) 2(1 - A) 

1 


X 


A(M2 - a!±!1 _ 
m(P, S')(/+m) ^ (/if + m)n(P, S') 


x2'Ki 


P- 


(AWfci) —ie (m^ +fcl) —ie 


(1—x)P+ 3^P+ 

When we perform the x integration, we have 


0/9 


$ 


= — I ag 6162 


X 


(A2+fc" , _ (Xl+(kx-r^r)- ie\ 

{l-x)P+ ' ' (a:-A)P+ 


(27) 


d?k 


(27r)2 P+ A (1 - A) 


X 


X 


A(M 2 - a!±P _ 

Il(P, S')(/+m) ^ (/if + m)M(P, S) 

1 

(m2 - (^a2 + (fcx 


/^r 


.2 (1-^) /• p(P‘^)(/+"^)J [{,}i + m)u{P, S) 

lag ^ 1^2 J ^2^^2 A2) 


where P is given in Eq. (^). 
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4.2.1 Calculation of/£^( A, rj^) 

One obtains f^{A, r^) from <f)“^ given in Eq. (|T^ by extracting the proton spin 
dependent part of 

4)"/^(7+)^“ = 2e^^Sirr{^ , (28) 


where = + 1 . 

We now apply this to 4*^ and for the numerator spinor contraction we calculate 


u{P, S){f+m)] \{^ + m)u{P, ^)1“ ( 7 +)^“ 


= Tr 


+ M)(- 75 ^)(/+ m) 7 +(/!f + m) 


= (__)(_4 ie+^^^) ^mP^^S^ik - ry + MPS^ik - 
= —2iP~^{AM + m)e'^^S^{kj_ — rj_y when x = A , 

i = —pLep^ypfja^lP&dP. 

I), we obtain 


(29) 


where we used = +1 and Tr 
When we insert Eq. (pOD into Eq. 


[S dependent part] 


(30) 


• 2 ( 1 -A) 

= — % ag 6162 


= — 2 ag^ 6162 

Then, from Eqs. 


r cPk± —2iP^{AM + m)e^^S^{k± — r±y 

P^{rl + B) J (27r)2 {kl + B) {{k^-r^y + XD 

(AM + m)(l - A) • f (fkj_ {k±-r^y 


{rl + B) 
and (^) we get 




{27ry{ki + B) {{k^-r^y + Xl) 


^ 1 2 M(AM + m)(l-A) 1 frl + B' 

7 r(A.rn - ^In 


B 


(31) 


From Eq. (|I^) we hnd that in terms of /i and the single-spin asymmetry 
transverse to the production plane in the SIDIS is given by (with St = y) 






_ n /it(A, rx) 2 

2 M/i(A,r^) M /i(A,r^) ^ 

Then, using the results in Eqs. (^) and (|^ , we get 


VyS^T = 


^ _ 6162 (AM-Fm)r^ r^ + B ^frf+B' 

/ 77 — : O 71 TT TTT O 


(32) 


Itt r]_ + {AM + my 
which agrees with Eq. (21) of Ref. p. 


B 


(33) 
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4.2.2 Calculation of/i^( A, rj^) 

Similarly, one obtains hi{A, r^) from given in Eq. ([T5|) by extracting the proton 
spin independent part of 

= 2r‘A , (34) 


where = |[7'^,7'^]. 

We again apply this to <hi? and for the numerator spinor contraction, we obtain 


^Y.HP<S){y+m)f[(li + m)u{P,S)]“ (<t7 


±s 

= -Tr 
2 


i+\l3a 


(Z’ + M)(/+ m)(-(7*7+ - 7’^7*))(/if + m) 


= —2i P~^ M(Afc® — xP) + m{k — r)* 

= —2i P~^{AM + m) {k± — r±y when x = A . 
Then, from Eqs. (^Tj), (0) and we obtain 


(35) 


,1.. s 1 9 M(AM + m)(l-A) 1 /r?+5\ 

^ - ag e,e, (36) 

Thus, from Eqs. (|3T|) and (^) we hnd the relation 

/{^(A,r^) = h^(A,r^). (37) 

is a special property of the quark-scalar diquark 

and (^^ schematically as 


We note that the equality Eq. 
model. 

We can write f^p and given in Eqs. 


/ir(A, ’’±) = ^^(A, r^) = (r^ + i?) ^^ ^ 


B 


with B as given in Eq. 


A = 


9 


and 
2 


2(27r) 


(- M (AM + m)(l - A) 

V 471- / 


(38) 


(39) 


We have the same formulas for and h-^ with A, r±, A, B replaced by A, r±, A, B. 

We note that we obtained /fZ and in Eq. (|3H ) from the hnal-state interaction 
diagram shown in Fig. |^(b). These are the functions relevant for semi-inclusive DIS 
0. The functions arising from initial-state interactions have an overall minus sign 
compared to those in Eq. (|38|), as pointed out by and conhrmed in [Q. However, 

f^rp and hi also have this property, therefore, the asymmetry factor u given in Eq. 

(pT|) is in fact independent of whether we use and from initial- or hnal-state 
interactions. 
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4.3 The cos 20 asymmetry 

We now consider the convolution terms in the numerator and denominator of the 
analyzing power v of the asymmetry (Eq. ([16D ): 


F = T 


(2h-pj^ h kj^ — 


G = T 


(fpj_(fk^S^{pj_ + kj_- q_L) (2 h pj^ h k^ - p^-k^ h^(A,p^)h^(A, k\_), 

fill 


= d p^d k^5 {p^ + k^ - g^)/i(A,p^)/i(A, fc^), 


(40) 


where we left out the flavor indices. With these dehnitions we can write 


V = 




■,2 p 

J- n 


M,M2Ea,aelGa 


(41) 


We will insert the schematic form Eq. (^2]) for fi and and Eq. (|3^ for hj- and hf. 
We hrst rewrite the denominator term G: 

G = y ^ exp {-^b^ ■ q^) MA, bl) j\{A, bl), (42) 

where we have dehned the Fourier transform of /i(A, k^) 

fi{A,bl) = y exp( 26 ^-p^) /i(A,p1) 

= 27rC' (^Koi^b) + , (43) 

where b = |6j^|, and similarly for /^. Thus, we obtain the exact expression for G: 

G = 27rGC r dbbJo(b\q,\) (Koi^/Bb) + ^^^^^bKl(^/Bb) 

Jo \ 2^/B 

X {Ko{^b) + . (44) 

Obtaining such an exact expression for F is much more difficult (if possible at 
all), hence we will express F in a form amenable to numerical evaluation. We hrst 
write 

F = -l -bMb\q^\)hi{A,b)hi{A,b), (45) 
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where we have dehned the Fourier transform 


hi{A,b) = [ exp {ih^ ■ p^_) h^(A,p^) 


= 2TiiA J dpJi{bp)\n 


'p‘^ + B' 


B Jp'^ + B' 

where p = |pj^|. This can be approximated from below by expanding 


In 


'p^ + B' 


p 


1 

+ X 


p 


1 

+ X 


p 


p"^ + B 1 2 + B I 3 + B 


+ ■■■ 


( 46 ) 


(47) 


For each term an exact Fourier transform expression can be obtained in terms of iF* 
functions. Keeping only the hrst term will lead for instance to 


h^{A,b) 


> 


1 'TT /4 ^ 1 r 

^ [KiiVBb) - -h^ [KoiVBb) + K2(VBb) 


(48) 


which is roughly a factor of 2 too small compared to the numerical evaluation without 
approximation. Eq. (^8]) leads to an asymmetry with approximately the right shape, 
but about a factor of 4 smaller in magnitude. This discrepancy can be reduced by 
taking further terms in the Taylor expansion into account. 

We will now investigate the obtained expressions for F and G by a numerical 
evaluation. In order to simplify the numerical calculation somewhat (since no absolute 
prediction can be made at this stage, because the overall magnitude of A and A are 
not known), we assume the situation of equal hadron masses (Mi = M 2 = M) and 
take momentum fractions such that B = B and D = D. This results in the following 
expressions, after expressing all dimensionful two-vectors in units of a/B, i.e. rescaling 
VBb^ and q_|_ — q^/y/B (idem for and kj_), 


2ttAA 

B2 J 

dbbJ2{b\q^\) 



2tiGC 

B J 

f dbbMb\q^\) 

(/'*pp'^”(''p)(p^+i)0 ■ 

(49) 


Next we make some generic choices for the various parameters. We take M = 
0.94GeV, m = 0.3 GeV, A = 0.8 GeV and D/B = 4, which implies that A 0.2 or 
0.5 and idem for A. 

Figure ^ displays the quantity P = BCCF/[AAG) as function of |gj^| in GeV 
(using A = 0.2, a/B Ri 0.24 GeV). The quantity P still has to be related to z/ which 
cannot be done without further assumptions. First of all, we will assume u quark 
dominance, which yields v ^ 2Fu/{MiM 2 Gu)- Next we will use some results obtained 
in Ref. 0, where the same asymmetry u was investigated and the following form was 
assumed (based on very general arguments and the simple model result of Ref. |Q): 


ftUA.pi) McMh 

/,(A,pi) ‘"pl + Ml;' 


(50) 
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Figure 6: Numerical result for P = BCCF/{AAG), using M = 0.94GeV,m = 
0.3 GeV, A = 0.8 GeV and A = A = 0.2. 

where Mh is the mass of the hadron and ch and Me were used as htting parame¬ 
ters. The values ch ^ ^ and Me ~ 2 were obtained from htting the 194 GeV data 
of the NAIO Gollaboration, by considering the case of one dominant quark havor 
contribution. In the present model calculation the ratio takes the form 

A l pI + B (p\ + B\ 

Unfortunately this shape is very different from Eq. (^OD for the choices of B and 
D made earlier (A 0.2 such that B k, 1/16 and D k, 1/4). Although both forms 
have similar large p\ behavior, it is mostly the small p\ behavior that is relevant. By 
comparing the curves resulting from Eq. (|50D (with c// = 1, Me = 2Mh = 2 GeV) and 
Eq. (^) (with D = AB = 1/4), one may expect 0.1 "^AjC <0.5, which then implies 
that 2 < |eie 2 | < 12 (incidentally this matches the value of |eie 2 | ~ 5, which was used 
for the numerical estimation in Ref. 0]). This range of values may then be used for 
crude estimates of asymmetries containing the function for a quark inside a proton, 
or equivalently an anti-quark inside an anti-proton. For a quark inside an anti-proton, 
or equivalently an anti-quark inside a proton, the overall prefactor is expected to be 
smaller. So if we restrict ourselves to pp collisions here and take A/C = A/C = 0.3, 
then one obtains as a very crude estimate: v ~ 2F/{M 1 M 2 G) = 2A‘^P/{BC‘^) ^ 3P, 
which means that the maximum of u is on the order of 30%. As said this is a very 
crude estimate and many assumptions went in. It cannot be viewed as more than an 
order of magnitude estimate, but we think it is an encouraging result. 

4.4 Discussion 

In this section we give a more general discussion of the qualitative features of the 
asymmetry, in particular its Q\ dependence. It may be good to note that the starting 
point of the calculation, that is, the factorized description of the asymmetry, requires 
that Q\ -C Q^, such that for large values of Q\ the asymmetry is not appropriately 
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described by the above formulas. At Q\ ~ Q^, the perturbative corrections will be 
the dominant source of an asymmetry. 

In order to obtain the general Q\ dependence of the asymmetry for small and 
large Q\, we start with the original convolution expression for F (the first line of Eq. 
(|i0|)). After multiplication by a trivial factor Q\IQ\ and using the integration to 
eliminate the delta function, we shift the integration variable p'^ = P± ~ 

to arrive at 


AA 


'Q 




J \2 


Q 




X- 


((p± + iQx? + B) - lqj_)2 + B 


(p'l + ili)^ + B ^ + 

In case B = B (which, as said, means equal masses of the initial hadrons and equal 
light-cone momentum fractions of the quark and antiquark), then one can perform 
symmetric integration to reduce 






+ p'lQ1 -2(gx 


Qi 

4 


(53) 


For small one can ignore the terms in the denominators and In terms of 

the expression Eq. (|5^), such that symmetric integration is appropriate and one can 
immediately conclude that F ~ Q\. 

Next we turn to the denominator of z/ (Eq. (p!6[) ) which is given by 


G 


CC J d^p'^ 


( p 1 + hQL? + D (pI - hQ±? + D 

iip± + + Bf (^(p'^ _ ig^)2 + 


(54) 


At Q\ = 0, G 7 ^ 0, hence the asymmetry [u F/G) vanishes. 

For large Q± the obtained model expressions do not yield an accurate description, 
although one does obtain a power law fall-off (see below), as one also would expect 
from perturbative QCD (which determines the large transverse momentum region). 
The point is that one runs into convergence problems. This applies for instance to 
the integral over fi (Eq. (pISD). Also does not fall off fast enough at large Q\ to 
guarantee convergence of certain Q^-weighted and integrated asymmetries (such as 
investigated in Refs, p, ^). Although one obtains a hnite result for 

J d^r^hi{A,rl) = (55) 

this is however not the object one encounters in the cos 20 asymmetry, nor in g^- 
weighted and integrated asymmetries. Rather one encounters in such weighted asym¬ 
metries the quantity 

J d^r±rlh^{A,rl), (56) 
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which diverges in the quark-scalar diquark model employed here. 

Therefore, one often assumes that the proton-quark-diquark coupling constants g 
are in fact form factors, see for instance Ref. [^. In the present quark-scalar diquark 


model calculation no such form factors are included (although the use of a regulator 
is implicitly assumed), because it would add another complication to the evaluation 
of the asymmetry and more importantly, in separating the perturbatively generated 
cos 20 asymmetry (which is only relevant at large Q'jJ, from the nonperturbative 

contribution x , one has to impose an upper cut-off on the Q\ range anyway. 
Our interest here is not in the specific fall-off of the asymmetry at large but 
rather in the moderate Q\ region, where the contribution to the asymmetry arising 
from initial-state interactions is maximal. 

For large Q\ one concludes from the above expressions (after including a regulator 
to insure convergence, e.g. a transverse momentum fall-off in g), that F,G and v ~ 
F/G decrease for large Q^. To see this in more detail, we will approximate F/G 
crudely by setting p_^ = 0 in the denominators and by ignoring B, B, D, D and the 
In terms altogether. In this way we obtain for the large Q± behavior of the ratio 


F 

G 


1 


(57) 


i.e. the asymmetry indeed falls off for large Q\. Since at small Q\ the ratio F/G 
grows as Q\, there has to be a turn-over in u as function of Q\, which has not yet 
been observed in experiment, but is clearly seen in the model calculation reported 
here. 

We want to emphasize that the quantities which determine the magnitude (and 
width) of the asymmetry u are the same as those appearing in the expression for the 
single-spin asymmetry proportional to hi x and in the context of the model also 
for the single-spin asymmetries discussed in Refs. |[ that depend on the Sivers 
distribution function. Thus, the parametric dependencies of these asymmetries can 
in principle be checked for consistency, in order to see whether it is at least consistent 
to assume that the asymmetries are generated by the same underlying mechanism. 
An example of such a comparison was given in Ref. . 

One final comment is on the scale. The model does not produce a dependence 
on that scale and the dependence of transverse momentum dependent asymmetries 
is a notoriously difficult problem (cf. e.g. |^, Q). Due to the lack of knowledge of this 

dependence, we can only expect that the asymmetry expression and the result from 
the model calculation should apply to the same range = (4 — 12 GeV)^) 

as that of the existing Drell-Yan data, from which we used fitting results. 


5 Conclusions 

In this paper we have studied the cos 20 distribution in unpolarized Drell-Yan lepton 
pair production within the context of a quark-scalar diquark model for the proton 
including an initial-state gluon interaction. Such initial- or final-state interactions 
lead to the appearance of (naive) T-odd distribution functions, such as the Sivers 
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effect function f^{x,p\) and its chiral-odd partner h^{x,p\) |^, We calculated 
those functions in the quark-scalar diquark model and found that they are equal in 
this model. Even though this equality is not expected to be satished generally in 
other models, this result does show that f^ixjp'jj and hj^{x,p\) are closely related 
and are expected to have similar magnitudes in general. With the model expressions 
for fi and we were able to write down an expression for the analyzing power u of 
the cos 20 asymmetry in the unpolarized Drell-Yan process. Under the assumption of 
u quark dominance and by using htting results of Ref. , we have given a numerical 
estimation of the asymmetry for the pp —> process. As an order of magnitude 

estimate we obtained for the maximum of u a. value of 30%. Despite the considerable 
uncertainty it is clear that based on this model calculation the cos 20 asymmetry 
can be of the same order of magnitude in pp —> as experimentally measured 

results in 7r“ Y —/i+ X (in the same range of values). It is natural to expect 
that the asymmetry in pp —> i~^i~X will be considerably smaller, but may still be 
expected to be on the percent level. 

Since the same mechanism (initial/hnal-state interactions) leads to nonzero func¬ 
tions f^{x,p\) and hj^{x,p\), it is clear that the single-spin asymmetries in the 
SIDIS and the Drell-Yan process are closely related to the cos 20 asymmetry of the 
unpolarized Drell-Yan process. Since the width and the magnitude of these asym¬ 
metries are determined by the same parameters in the model, one can relate the 
asymmetries and this may be tested by experimental data. All this provides new in¬ 
sight into the role of quark and gluon orbital angular momentum as well as of initial- 
and hnal-state gluon exchange interactions in hard QCD processes. 
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Appendix: Derivation of Eq. (|23) 

We present the derivation (based on Ref. |3) of the starting formula of Eq. 
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u{q + r)Y 


^ ^ — M 
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{k — r)+ -|- ie 
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— u{q + r )7 {k — r)’ 


— u{q + r) 2 


^ ^ — M 
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'+d-M 
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{k — r)+ -f- ie 

— - ^ . u{q + r)[{^+ 4-M)-{4+f-M) 

[k — ry + le l 

2 1 
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- 4 ~\~ 4 — ^ 
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{k — r)+ -f ie 


+ 4-M 
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(58) 


2 

{k — r)+ + ie 


u{q + r)'j^ 


5 


where we used the equation of motion u{q + r){^+ / — M) = 0 in the fourth line. In 
the above, 2/((fc — r)~^ + ie) is the eikonal propagator. 

Before going from the second to the third line in Eq. ([58|), we deformed the 
contour of integration to the upper inhnity in the complex {k — r)~^ plane so that 
\{k — r)+| S> \{k — 'r)“|, \{k — r)*| is satisfied along the new contour. We note that 
what we deformed is the line along which the {k — r)~^ integration is performed, and 
the pole position of {k — r)^ + ie = 0 (at which we compute the value of the residue) 
is not influenced by this deformation. 


References 

[1] D. Sivers, Phys. Rev. D 41, 83 (1990); Phys. Rev. D 43, 261 (1991). 

[2] Z. Liang and T. Meng, Phys. Rev. D 42, 2380 (1990); Phys. Rev. D 49, 3759 
(1994); 

C. Boros, Z. Liang and T. Meng, Phys. Rev. Lett. 70, 1751 (1993). 

[3] J. Qiu and G. Sterman, Phys. Rev. Lett. 67, 2264 (1991); Nucl. Phys. B 378, 52 
(1992). 

[4] J.C. Collins, Nucl. Phys. B 396, 161 (1993). 

[5] M. Anselmino, M. Boglione and F. Murgia, Phys. Lett. B 362, 164 (1995). 

[6] D. Boer, Phys. Rev. D 60, 014012 (1999). 

[7] S.J. Brodsky, D.S. Hwang and 1. Schmidt, Phys. Lett. B 530, 99 (2002). 

[8] S.J. Brodsky, D.S. Hwang and 1. Schmidt, Nucl. Phys. B 642, 344 (2002). 

[9] J.C. Collins, Phys. Lett. B 536, 43 (2002). 

[10] D. Boer, |hep-ph/0206235| . 

[11] D. Boer and P.J. Mulders, Phys. Rev. D 57, 5780 (1998). 

[12] G.R. Goldstein and L. Gamberg, |hep-ph/0209085| . 

[13] NAIO Collaboration, S. Falciano et ai, Z. Phys. C 31, 513 (1986); 

NAIO Collaboration, M. Guanziroli et al, Z. Phys. C 37, 545 (1988). 

[14] J.S. Conway et a/., Phys. Rev. D 39, 92 (1989). 

[15] K. Hagiwara, K. Hikasa and N. Kai, Phys. Rev. D 27, 84 (1983). 

[16] T. Gehrmann, |hep-ph/9608469| ; 

M. Ahmed and T. Gehrmann, Phys. Lett. B 465, 297 (1999). 


21 










[17] C.S. Lam and W.K. Tung, Phys. Rev. D 21, 2712 (1980). 

[18] A. Brandenburg, O. Nachtmann and E. Mirkes, Z. Phys. C 60, 697 (1993). 

[19] J.C. Collins and D.E. Soper, Phys. Rev. D 16, 2219 (1977). 

[20] A. Brandenburg, S.J. Brodsky, V.V. Khoze and D. Muller, Phys. Rev. Lett. 73, 
939 (1994). 

[21] K.J. Eskola, P. Hoyer, M. Vantinnen and R. Vogt, Phys. Lett. B 333, 526 (1994); 
J.G. Heinrich et ai, Phys. Rev. D 44, 1909 (1991). 

[22] E.L. Berger and S.J. Brodsky, Phys. Rev. Lett. 42, 940 (1979); 

E.L. Berger, Z. Phys. C 4 (1980) 289; Phys. Lett. B 89, 241 (1980). 

[23] J.P. Ralston and D.E. Soper, Nucl. Phys. B 152, 109 (1979). 

[24] R.J. Fries et ai, Phys. Rev. Lett. 83, 4261 (1999); Nucl. Phys. B 582, 537 (2000). 

[25] D. Boer, P.J. Mulders and O.V. Teryaev, Phys. Rev. D 57, 3057 (1998). 

[26] A.V. Belitsky, X. Ji and F. Yuan, |hep-ph/0208038| . 

[27] G.T. Bodwin, S.J. Brodsky and G.P. Lepage, Phys. Rev. D 39, 3287 (1989). 

[28] R.L. Jaffe and X. Ji, Nucl. Phys. B 375, 527 (1992). 

[29] R. Jakob, P.J. Mulders and J. Rodrigues, Nucl. Phys. A 626, 937 (1997). 

[30] S.J. Brodsky, D.S. Hwang, B.Q. Ma and 1. Schmidt, Nucl. Phys. B 593, 311 

( 2001 ). 

[31] X. Ji and F. Yuan, Phys. Lett. B 543, 66 (2002). 

[32] D. Boer, Nucl. Phys. B (Proc. Suppl.) 79, 638 (1999). 

[33] D. Boer, Nucl. Phys. B 603, 195 (2001). 

[34] A.A. Henneman, D. Boer and P.J. Mulders, Nucl. Phys. B 620, 331 (2002). 


22 




